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Abstract 

Semiclassical quantization is exact only for the so called solvable potentials, such as the harmonic 
oscillator. In the nonsolvable case the semiclassical phase, given by a series in h~, yields more or 
less approximate results and eventually diverges due to the asymptotic nature of the expansion. 
A quantum phase is derived to bypass these shortcomings. It achieves exact quantization of 
nonsolvable potentials and allows to obtain the quantum wavefunction while locally approaching the 
best pre-divergent semiclassical expansion. An iterative procedure allowing to implement practical 
calculations with a modest computational cost is also given. The theory is illustrated on two 
examples for which the limitations of the semiclassical approach were recently highlighted: cold 
atomic collisions and anharmonic oscillators in the nonperturbative regime. 
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The semiclassical treatment of integrable systems, which can be traced back to Bohr's 
atomic model of planetary motion and Einstein's classic paper on the quantization of reg- 
ular motion 1] is assumed to be a well-established and venerable subject. It is true that 
WKB theory, where the phase to be quantized is the classical action, can be found in any 
standard textbook, but WKB often results in approximations that are quantitatively too 
crude and that fail to capture the physics of the problem. Indeed, WKB theory achieves 
exact quantization for a restricted number of potentials - known as the solvable potentials 
- such as the harmonic and Morse oscillators or the centrifugal Coulomb problem. In the 
last decade, the application of supersymmetric (SUSY) methods to quantum mechanics has 
enlarged this list to a handful of other potentials, quantized by employing SUSY WKB Q|. 
However, even in the solvable cases, the WKB wavefunctions are innaccurate especially at 
the turning points where they blow up, and a consistent divergence-free WKB scheme is still 
a topic of current investigation g. In the more general nonsolvable case, WKB quantiza- 
tion has frequently resulted in useful approximations to compute the energy levels of excited 
states, but recently several shortcomings were pointed out: for example the phase loss in the 
classically forbidden regions is badly taken into account by WKB theory in potentials used 
in atomic clusters calculations [4j; in cold atom collisions WKB quantization breaks down 
for very excited states |5[; for anharmonic potentials, the failure of WKB has prompted 
extensive developments of numerically involved quantum techniques with the aim of obtain- 
ing accurate results [f], Q for applications ranging from molecular physics to quantum field 
theories. 

We show in this work that these shortcomings can be resolved by considering an ex- 
act quantum phase. Although the quantum phase bears a very close relationship to the 
semiclassical one, it is necessary to go beyond the semiclassical approach in the presence of 
nonsolvable potentials. Indeed, the failure of the semiclassical approximation resides in the 
limitations imposed by the underlying classical dynamics. Let us take a one dimensional 
conservative system at energy E, which will be our main concern in this work. The WKB 
wavefunction is built from blocks given by s| 

V>wkb(x, E) = A{x, E) exp[iS(x, E)/h] (1) 

where the phase function S(x,E) is the classical action and A 2 (x, E) = (c^S*)" 1 gives the 
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classical probability amplitude. The WKB quantization condition reads 

S{t 2 ,E)-S(t 1 ,E) = (n+^7th, (2) 

where tip are the turning points, \i is the Maslov index (generally 2 in one dimensional 
systems) and n is the level integer. To improve the approximation a semiclassical expansion 
going beyond these purely classical terms can be carried out by going to higher order in h. 
However this expansion is an asymptotic series, which means that although going to higher 
order may improve the accuracy of the results, at some point the series generally diverge 
0| . Moreover the divergence of the amplitude at the turning points becomes worse at each 
order and regularization techniques are needed to compute the phase integrals. 

This is why a quantum phase a(x, E) while implicitly summing the divergent semiclassical 
expansion needs to be defined from the start by appropriately transforming the Schrodinger 
equation 

h 2 d 2 x y(x)+p 2 (x)y(x) = J (3) 

where p(x) is the classical momentum. A transformation of the Liouville- Green type is 
taken by writing 

y(x) = (d x ax)r 1/2 w(0, (4) 

so that appears as a 'phase' and the prefactor as an amplitude, as in Eq. (0). Assume 
w(£) fulfills the equation 

h 2 6lw(Z) + R(Z)w(Z) = 0, (5) 
where the choice of the unspecified function R(£) determines the choice of w. £ then obeys 

fl(0(fl , »e) 2 -p a (^) + y<e;*> = o, (6) 

where (£;x) = d x £,/d x £, — ^{d 2 ^/d x ^) 2 denotes the Schwartzian derivative. We thus see that 
the choice of the phase first depends on the choice of the carrier function w (or equivalently, 
of -R(O)) an d then on the choice of the boundary conditions that need to be imposed on 
the third order nonlinear Eq. (JBj). This is an illustration of the ambiguity suffered by phase 
functions in quantum mechanics, due here to the fact that there is no unique manner to cut 
a given wavefunction into a phase £ on the one hand, and an amplitude function obeying 
the continuity equation a(x,E) = {d x ^)~ 1 / 2 (as in the semiclassical case) on the other. 

Formally can be expanded as an asymptotic series in h irrespective of the specific 
choice of R. It is nevertheless apparent from Eq. (JHJ) that the only choice that will give 



£(x) ~ S(x) to first order in H corresponds to R(£) = ±1, leading by Eq. (JSJ) to circular 
or exponential carrier functions. The price to pay is that the H expansion based on these 
functions - such as the WKB approximation, necessarily diverges at the turning points. 
However the exact solution of Eq. (JHJ) does not present such deficiencies, pointing to the 
possibility of defining a quantum phase from Eq. (0) with R(£) = ±1. For bound state 
problems the solutions of Eq. (J3J) are real so w can be taken as a circular function. y(x) 
then becomes proportional to ct(x) sina(x), with o denoting the phase £ when R = ±1. 
Let p(x) span the interval ]si,s 2 [ (typically ] — oo,+oo[ or ]0, +oo[ for radial problems). 
It can then be shown that a is a positive definite quadratic form that behaves as 
a(x — > Sj = i j2 ) — > oo so that setting cr(si) = the quantization condition reads 

a(s 2 ,E) = (n + l)n, (7) 

where n is the level integer as in Eq. Eq. (jJJ) is exact and holds irrespective of the 

boundary conditions imposed on a . 

For computational purposes, Eq. (JBJ) is intricate to solve. It can be checked that by 
writing (we use atomic units and assume R(£) = 1 from now on) 



M{x, E) = d x 



1 

a(x,E) + -ln(d x a) 



(8) 



Eq. (JBJ) leads to the complex but first order nonlinear differential equation 

d x M = i (p 2 (x) - M 2 (x)) = F(M{x), x), (9) 

where T denotes the middle term taken as a functional. Our strategy will consist in solving 
the equation for M; the real part will give us d x a, which can be numerically integrated to 
obtain a, yielding both the wavefunction and the total phase a(s2,E). To do so we first 
linearize the equation for M by expanding the functional to first order in the vicinity of an 
initial trial function Mq(x). We then solve 



d x M q+1 = F(M q (x),x) + 
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(M q+1 (x) - M q (x)) , (10) 



with q = and where 5 stands for the functional derivative. Eq. is a linear first 

order differential equation that can be solved straighforwardly. Of course since T has been 
linearized Eq. (fTUj) is not equivalent to Eq. (JSJ), hence the subscript in M\. Convergence 
towards M is achieved by iterating the procedure, i.e. we now solve Eq. (fTUj) for q = 1, 



obtaining a better approximation M2 and so on. This iterative linearization procedure, 
known as the quasilinearization method (QLM) replaces the solution of a nonlinear differen- 
tial equation by iteratively solving a linear one. It was introduced 3 decades ago by Bellman 
and Kalaba |l2| in the context of linear programming, but the extension of QLM to the 
type of functions dealt with in quantum mechanics is quite recent [13]. In particular, the 
property of quadratic convergence which makes this method powerful still holds. Indeed, 5 
to 6 iterations typically suffice to obtain solutions with high precision (e.g. more than 20 
decimal digits). 

To iteratively solve Eq. ([ID}) , two ingredients are needed: first the trial function Mo(x), 
second, the boundary condition, M q (xb) which is set from the start since it must be the same 
for each q. The choice of M is actually unimportant, since any sufficiently well behaved trial 
function, i.e. behaving as \p(x)\ between the turning points and as i \p(x)\ (resp. — i \p(x)\) 
for x < ti (resp. x > t 2 ), with a smooth transition between the regions, will lead to the same 
converged answer. An elegant solution consists in determining the first order semiclassical 
expression not diverging at the turning points, which is done by appropriately choosing R(£ ) 
in Eq. (jSJ) and expressing a in terms of but we shall not pursue this task here jl4j |. 

What is crucial is the choice of the boundary condition because it determines the behavior of 
the solution. We draw here on previous work where we had shown that if a(x) is written as 
an infinite h expansion, there is a unique boundary condition that yields the classical action 
S(x) when the limit H — > is taken For other boundary conditions, the limit h —> 
yields a semiclassical phase (and the corresponding amplitude) that displays the oscillations 
of the WKB wavefunction. The existence of an optimal boundary condition is well known in 
the context of solvable potentials, where the use of analytic solutions (e.g. parabolic cylinder 
unctions for the harmonic oscillator, Whittaker functions in the centrifugal Coulomb case 
I5I ]) allows to explicitly construct a nonoscillating quantum phase. However for nonsolvable 
potentials a pro cedure based on special functions does not exist, whereas the h expansions 
employed in [ll| are of formal nature. To ensure that the quantum phase behaves as the h 
expansion would if it converged we choose here the simplest solution, namely to pick 
an appropriate point x^ in the classically allowed region where the standard (i?(£) = 1) 
semiclassical expansion can be employed and go to the highest possible order in h before the 
series starts to diverge. More refined methods of estimating the optimal boundary condition, 
based on super and hyper-asymptotic expansions are in principle available, but they 
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FIG. 1: Left scale: The quantum phase cr(x) (solid black curve) for the anharmonic oscillator at 
E = 5 (units in a.u.) shown along with the classical action S(x) (thick gray/red curve) plotted 
between the turning points (semicircles on the x axis). Note that we have taken S(ti) = 7r/4. 

n 

This 7r/4 constant is familiar from standard WKB analysis |8( but it can also be rigorously proved 
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that the semiclassical phase not diverging at the turning points, obtained from Eq. © 
with i?(£) = £, behaves as cr sc i(x) ^ S(x)/h + 7r/4 for x ^> t\ (but o sc i(t\) = vr/6). Right scale: 
The associated wavefunction asino". Since E is not an eigenvalue this improper wavefunction is 
normalized per unit energy, as usual in scattering theory. 

are unnecessary insofar as neither the behavior of a nor the quantum mechanical quantities 
would be affected by their use. 

We will now illustrate and detail the properties of the quantum phase on two specific 
examples. The first example concerns anharmonic oscillators which are employed in the 
investigation of quite different phenomena, ranging from molecular vibrations to quantum 
field theories and phase transitions. This has generated different schemes to compute the 
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or delicate resum- 



eigenvalues and the wavefunctions, calling for large numerical basis 
mation techniques The failure of the semiclassical approximation which is important 
for the lowest levels (e.g., for symmetric potentials x 2m WKB quantization gives the wrong 
behavior of the energies as a function of m) but persists for more excited states, has been 

n 

generally attributed to the existence of complex turning points away from the real line [19]. 
However the complex plane does not play directly any role in the present scheme, based on 
the construction of a real quantum phase. Fig. 1 shows the quantum phase a(x, E) along 
with the classical action for a Hamiltonian with an anharmonic potential x 2 + 2x 4 . Given the 
symmetry of the potential, we have taken Xb = and fixed the value Mo(xb) by carrying out 
a semiclassical expansion up to o(h 14 ). We have also plotted on the same figure the associ- 
ated wavefunction. It diverges beyond the right turning point because we have voluntarily 
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FIG. 2: Energy levels of the anharmonic oscillator (a.u.) obtained by quantizing the quantum 
phase. The smooth curve represents a(s2,E)/ir as a function of the energy and the boxes are 
plotted each time this quantity equals the integer n + 1, determining the exact levels (here the first 
20 levels n = — 19 are plotted). The inset zooms on the ground state: the thick gray/red line 
gives the semiclassical quantization curve (S(t2,E) — S(t\, E))/tt + 1/2 (the H expansion diverges 
beyond this order) and the triangle the semiclassically quantized level, which is off by more than 
10%. 

chosen E not to be an eigenvalue (here E lies between the 2nd and the 3rd levels). The 
bound states are found by applying the quantization condition (J7J) in the following way: the 
total phase a(s2,E) at s 2 = oo is determined for different values of the energy (the energy 
grid can be more or less sparse, depending on the required numerical precision). The points 
are then interpolated to yield the curve a(s2,E) as a function of E seen on Fig. 2, which 
allows to solve for E in Eq. (JJJ). By that method we have determined the ground and the 
first 35 states of this anharmonic oscillator and compared the eigenenergies to the precise 
quantum calculations of obtaining exactly the same results as those published in their 
Table III. 

The second illustration involves a potential well with a strong repulsion at short-range and 
a long-range attractive tail. This type of potential is of interest in the study of cold atomic 
collisions, a field that has been sparked by the development of photoassociative spectroscopy 
0|. The WKB approximation breaks down for excited states near the threshhold, a fact that 
is hardly surprising since it is understood that WKB may fail when the quantum particle 
explores large areas of classically forbidden regions j^J. We take the following 12-6 Lennard- 
Jones (LJ) potential with the classical momentum given in the scaled form as 



p 2 {x) = B 
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FIG. 3: Derivative of the quantum phase d x a (solid black curve) and of the classical action d x S = p 
(thick gray/red curve) in the 12-6 LJ potential characteristic of cold atom collisions (scaled a. u.). 
These quantities, plotted for the most excited bound state, are shown in the zone near the outer 
turning point ti ~ 9.5. The inset shows the zone near the inner turning point t\ ~ 0.9 (note the 
very different scales). 

where B is a 'strength' parameter. This potential with B = 10 4 has often been employed as 



a benchmark (see 
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21 1 and Refs. therein) and is known to support 24 states. The WKB 
quantization condition (J2J) gives energies with an error relative to the local level spacing that 
globally increases with E. Fig. 3 compares the derivative of the quantum phase d x a(x, E) 
with the classical momentum d x S(x, E) for the last bound state just below the threshold. Xf, 
was taken to the right of the potential minimum and the h expansion was carried out up to 
the 12th order. The 2 curves are barely distinguishable in the classically allowed region, so 
we have focused on the zone near the turning points. Note that the quantum curve largely 
penetrates into the classically forbidden region well beyond the outer turning point. This is 
typical of excited states in potentials with a long-range attractive tail and can be seen as the 
underlying reason ruling the breakdown of the semiclassical approximation. Quantization 
proceeds as in the first example, by determining a(s2, E) on an energy grid then solving for 
Eq. (JZj). The resulting curve is shown in Fig. 4 and the quantized energies of the 24 levels 



exactly reproduce the values of exact quantum mechanical calculations [20| . In particular 
the position of the last level determines the scattering length, a crucial parameter in the 
production of Bose-Einstein condensates. 

To summarize, we have employed a quantum phase that goes beyond the semiclassical 
approximation in giving the exact quantization energies as well as the wavefunctions. We 
have also given a numerical procedure to determine the phase based on the linearization of 
the phase equation and a boundary condition obtained from a local semiclassical expansion, 
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FIG. 4: Exact energy levels (boxes) of the 12-6 LJ potential obtained by quantizing o~(s2,E) (solid 
curve), in scaled a.u. The inset zooms on the last state below threshold (E = 0) and also shows 
the WKB prediction, about 40% too low (triangle). 

and illustrated the approach in the case of two nonsolvable potentials. 
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